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(N| Abstract 



The effect of magnetic impurities on the ballistic conductance of nanocontacts is, as suggested in recent work, amenable to ab initio 
^SJ study |Q]. Our method proceeds via a conventional density functional calculation of spin and symmetry dependent electron scat- 
tering phase shifts, followed by the subsequent numerical renormalization group solution of Anderson models - whose ingredients 
and parameters are chosen so as to reproduce these phase shifts. We apply this method to investigate the Kondo zero bias anomalies 
that would be caused in the ballistic conductance of perfect metallic (4,4) and (8,8) single wall carbon nanotubes, ideally connected 
to leads at the two ends, by externally adsorbed Co and Fe adatoms. The different spin and electronic structure of these impurities 
1—1 are predicted to lead to a variety of Kondo temperatures, generally well below 10 K, and to interference between channels leading 
_ to Fano-like conductance minima at zero bias. 
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1. Introduction 

The zero-bias anomalies observed in STS conductance spec- 
troscopy through adsorbed magnetic impurities and to some ex- 
tent in metal break junctions have recently revived interest in 
the Kondo effect. Addressing these systems theoretically poses 
several problems. In the first place, and unlike quantum dots, 
ab initio electronic structure calculations such as density func- 
tional theory (DFT) are essential to establish a quantitatively 
meaningful starting point. Which among the impurity-related 
levels and resonances drive the spin polarization, what is their 
multiplicity, their hybridization, etc. are all questions that need 
an ab initio calculation. Next, this information must be trans- 
lated into some manageable many body Hamiltonian, possibly 
without the loss of the brute quantitative information provided 
by DFT. Finally, the many body Hamiltonian(s) must be solved, 
to extract Kondo parameters and the predicted conductance fea- 
tures near zero bias, possibly with their behavior with param- 
eters such as nanocontact geometry, temperature and external 
field, to be eventually compared with experiment. One ap- 
proach in this direction was recently taken by our group Q. 
Given a nanocontact between two leads, one identifies, with the 
help of symmetry, the impinging and outgoing channels that 
carry current across the impurity. From the matching symme- 
try selected local densities of states at the impurity, one identi- 
fies the important impurity orbitals with their different magnetic 
splittings and hybridizations. This leads to formulate multi- 
orbital Anderson models, which contain a multiplicity of pa- 
rameters to be adjusted. In our scheme the parameters are 
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adjusted to yield, within the Hartree-Fock approximation, the 
same channel- and spin-dependent impurity scattering phase 
shifts as those that we calculate ab initio by DFT - whose input 
information is therefore put to maximal use. For the last step, 
solving the Anderson models, we employed a standard numer- 
ical renormalization group (NRG) scheme. While other groups 
have dealt with the overall Kondo problem in different ways 
\2, 3 1, we find our "DFT + NRG" route extremely instructive, 
and worth exploring in more complex situations than the simple 
Au-Ni-Au contact studied in Ref (T|. In the present application 
we consider a single wall carbon nanotube (SWNT) as our lin- 
ear conducting system, and a single externally adsorbed transi- 
tion metal atom, either Co or Fe, as the magnetic impurity. To 
begin with, the metallic nanotube has two conducting channels 
instead of only one as Au. The magnetic atoms in turn have in 
principle a richer multiplicity of magnetic levels than Ni. We 
wish to explore what this richness might bring. 

2. Systems and symmetries 

We considered alternatively Co or Fe impurities on either 
(4,4) or (8,8) metallic SWNTs (see fig. [T}. If z is the SWNT 
axis, its electronic states of can be classified according to parity 
with respect to xy plane reflection (e-o, even- odd) and xz plane 
reflection (s - a, symmetric and antisymmetric). DFT calcula- 
tions (see section [3]l predict that the externally adsorbed impu- 
rities should have minimum energy when at the hollow site (see 
fig. |2j, that is above the center of a carbon hexagon. Assum- 
ing that geometry, the impurity electronic states can be clas- 
sified according to the same parity numbers as those of clean 
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SWNTs. We are interested in particular in 3d and As impurity 
orbitals, whose parities are shown in tab. [T] 

3. Ab initio electronic structure 

We carried out standard density-functional theory (DFT) cal- 
culations, allowing for full relaxation of all atomic positions 
in a unit cell, which comprised 80 and 160 carbon atoms for 
the (4,4) and (8,8) tubes respectively plus one Co or Fe ad- 
sorbed impurity. Calculations used the standard plane-wave 
package Quantum-ESPRESSO H within the generalized gra- 
dient approximation (GGA) to exchange-correlation function- 
als in parametrization of Perdew, Burke and Ernzerhof. The 
plane wave cut-offs were 30 Ry and 300 Ry for the wave func- 
tions and for the charge density, respectively. Integration over 
the one-dimensional Brillouin zone was accomplished using 8 
k-points and a smearing parameter of 10 mRy. When necessary 
to test the sensitivity of DFT results to correlation effects, we 
extended to "GGA+U" with a reasonably small Hubbard "U" 
[4|- but generally the straight DFT result was used. 

We found that Co behaves as a S = I impurity on both (4,4) 
and (8,8) SWNTs, its d xz orbital driving the spin polarization. 
The Co atom switches from the 3d 1 As 1 configuration of the iso- 
lated atom to a slightly surprising low-spin 3d 9 4s° one when 
adsorbed on the nanotube. Fe behaves as a S = 1 impurity on 
the (8,8) tube, similarly switching from the high-spin 3d 6 4s 2 
of the isolated atom to a low-spin 3d & 4s° in the adsorbed state. 
Here the pair of orbitals d xz and d xv is magnetically polarized 
(see fig. [3]and tab. [2]), a result in good agreement with previous 
calculations ||5][6). Orbitals d xy for Co and d z i for Fe are partly 
empty, and fall near the Fermi energy in straight DFT: but they 
promptly move below Ef when even a small U is switched on. 
We conclude that these orbitals are not going to be involved in 
Kondo behaviour and can be neglected to a first approximation 
in order to keep the many-body model simple. The behavior of 
Fe/(4,4) is complicated. The s orbital is partly filled, and d z i 
is magnetically polarized besides the (d xz ,d xy ) pair, so here Fe 
should behave as a S — | impurity. 

As in previous work U we implemented DFT computa- 
tion of the (spin-polarized) mean-field ballistc conductance and, 
more importantly, of the impurity-related spin- and channel- 
selected phase shifts suffered by the SWNT conduction elec- 
trons as a function of energy. An example is shown in fig. |4]for 
Co on the (4,4) SWNT. 



4. Generalized Anderson model 

The Kondo model is usually understood by means of a many- 
body Anderson Hamiltonian [7|. In our case we need to extend 
it in principle to the four SWNT conduction bands, each hy- 
bridized with some impurity orbital among the 3d and 4s, of 
same symmetry. These impurity orbitals in turn are mutually 
coupled by an intra-atomic ferromagnetic Hund exchange term, 



H = 
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(3) 
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where d\ creates an electron on the impurity orbital with sym- 
metry i, ct creates an electron in a k conduction state with sym- 
metry i, sare the Pauli matrices, f, is a potential scattering term 
due to to the charge density of the impurity, V,- is the coupling 
of the impurity orbital with the conduction electron states, e, is 
the bare energy of the impurity orbital, [/, is the Hubbard repul- 
sion on the orbital, J < is a global Hund exchange parameter 
(favouring high spin for the isolated impurity), and the single 
particle energies of conduction electrons e# are such as to give 
a constant density of states, with exactly the same value as that 
of the clean SWNT (per spin direction) as computed by DFT: 

6{e — — p — - — — , i = es, eo,as, ao (5) 
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In practice, only conduction bands coupled to a magnetic or- 
bital are retained in our NRG procedure (see section [SI. This 
leaves us with a single band coupled to a single impurity level in 
the case of Co (once orbital d xy is ignored), and with two bands, 
each coupled to one impurity level, in the case of Fe/(8,8) (once 
orbital d z i is ignored). The case of Fe/(4,4) is more involved and 
we will presently not deal with it. 

5. Joining up DFT and many body 

Hamiltonian eq. [T] is easily solved in the (unrestricted) 
Hartree-Fock approximation [7], breaking spin rotational sym- 
metry. This leads to a phase shift in conduction electrons of 
symmetry i (i = es, eo, as, ao) at the Fermi energy 

F, 



df = (pi + arctan ■ 



(6) 



where <pi = arctan npti =a is the phase shift caused by the im- 
purity charge scattering. This is numerically found to be neg- 
ligible, so we shall ignore it from now on. The peak of the 
impurity DOS is found to be at 



e'[ = E, + U i (n i ,T ) - o--(m tot - m,) 
where 

1 T a 
(nf) = - arctan — 
7i e\ 

is the average occupation of up/down orbital 

rrii = (nj) - {n\) 
is the magnetization of each orbital and 
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is the total magnetization of the atom. As in |Qj, we choose to 
reproduce the same phase shifts at the Fermi energy for each 
symmetry, and the same peaks in the density of states of the 
impurity orbitals as those computed by DFT. This allows to 
uniquely fix e,, t/, and F, as long as just one magnetic orbital 
is considered in eq. [T|- that is the case of Co (i = os). When 
more than one orbital is involved, such as in Fe, (or in Co if or- 
bital d xy were to be taken into account) we need to fix J as well. 
We can extract J from the DFT calculated exchange splitting of 
filled orbitals, according to 
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(11) 



Since different d orbitals have slightly different splittings, we 
just took an average value as deduced form different orbitals. 
This yields J ~ 1 eV in Co and J ~ 1.2 eV in Fe. 

6. Results of NRG calculations 

We solved the Anderson Hamiltonian by means of NRG 
1 8 9 1, which allows to compute all the needed static and dy- 
namic quantities we need in an almost exact, albeit numeri- 
cal, way. We extracted the conduction electrons phase shifts 
from the single particle energies at the zero energy fixed point, 
and the Kondo temperature from the impurity Green function at 
imaginary frequency: 
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(Z par , is the partition function and GS the ground state). The 
Kondo temperature is given by 



nwZr 
4k b 



(13) 



where w = 0.4128 is the Wilson coefficient and Z is the quasi- 
particle residue 



2-1 = 1 dW 
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(14) 



Alternatively, an approximate formula [10], valid for one impu- 
rity coupled to one channel, 



0.4107 



UiT, 



(15) 



could be used, with similar results. 

The zero-bias conductance is given, in terms of the final 
phase shifts, by 

g=3~= cos 2 ((5 f J - 5 0S ) + cos 2 (6 ea - 5 oa ) = g s + g a (16) 
(jo 

where Go = =jr is the quantum of conductance. Note that in 
the clean tube G = 2Gq. Phase shifts are only computed for 
Kondo channels, and are found to be always ^ jt/2. For the 
non-Kondo channels, they can be directly extracted from DFT. 



Since in DFT they are ^ 0, they can be safely neglected. Sum- 
ming up, both Co/(4,4) and Co/(8,8) should exhibit a (zero tem- 
perature and zero bias) conductance G ~ Go, whereas Fe/(8,8) 
should have G ~ 0. These results remain valid so long as ei- 
ther temperature and/or bias remain well below Tk- However, 
it turns out that Kondo temperatures 7j are quite low (see tab. 
[3J, which might make this effect hard to observe in a real ex- 
periment. Interestingly, a much higher Kondo temperature of 
about 15 K has been quoted for Co/ graphene ifTTI . While the 
reasons for this difference between graphene and nanotubes are 
presently being investigated, it should be noted that several fac- 
tors differ, including symmetry, and heavy doping in real, de- 
posited graphene. 

Finally, we can qualitatively address the predicted bias- 
dependent lineshape of the Kondo conductance anomaly. 
Through the Keldysh technique for non-equilibrium Green- 
functions it is possible to compute the finite-bias conductance 
Ifl2l . once the impurity Green function G,(e) is calculated from 
NRG El: 



g s , a = 1 - r3G,(e) 
For simplicity, we have taken 



G,(e) 



ryr 

e + iT k 



where 
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k b T K = ~^T K = 0.342F, 
This gives rise to a Fano lineshape IfMl 



(q + vf 



( 9 2 + l)(v 2 + l)' 



v = e/T k 



(17) 



(18) 



(19) 



(20) 



with q — 0, so for each band (s - a) the lineshape is predicted 
to be a symmetric antilorentzian, with a width proportional to 
the Kondo temperature. Small sources of asymmetry will arise 
from a) the potential scattering f, which we ignored in Eq. |4jb) 
from the interference with orbitals belonging to the same band 
a or s, but with different symmetry e/o; and c) from particle- 
hole asymmetries in eq. 18 However, we estimate that the 
asymmetry parameter q should generally remain below 0.1. In 
Co/(4,4) and Co/(8,8), only g s contributes to the lineshape, g a 
being almost one - and moreover independent from energy on 
the Kondo energy scale. In Fe/(8,8), both g s and g a have an 
antilorentzian shape, although with very different widths. The 
total lineshape is just their sum (see fig. [5j. 

7. Conclusions 

We implemented our recently devised DFT+NRG scheme 
CI to calculate the Kondo effect caused by Co and Fe adsorbed 
impurities on the conductance of (4,4) and (8,8) nanotubes. On 
the methodological side, the present calculation represents a 
good pedagogical illustration of our technique. For the systems 
chosen, the predicted anomalies are symmetric antilorentzian 
dips, reducing total zero bias conductance to zero for Fe, and 
by a factor 1/2 for Co. While there are no data to compare with, 
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Table 1 : Symmetries of d and s orbitals with respect to the xy plane (even e - 
odd o) and the xz plane (symmetric * - antisymmetric a). 



Impurity 


Magnetic Orbital 


Symmetry 


Spin 


Co 


dxz 


OS 


l 

2 




(dxy) 


ea 





Fe 


dxz 


OS 


1 

2 




d xy 


ea 


1 

2 




(d z 2,s) 


es 






Table 2: Magnetic orbitals as found from DFT calculations, their symmetry, and 
spin they carry (S = 1/2 for each magnetic orbital). Orbitals in parentheses are 
not Kondo orbitals, so do not contribute to the total spin of the impurity and are 
ignored in the many-body model, but still participates in transport. 



Imp. 


SWNT 


Orb. 


S/(eV) 


U(eV) 


r(eV) 


7Y(K) 


Co 


(4,4) 


dxz 


-1.62 


2.17 


0.082 


~ 1 




(8,8) 


dxz 


-1.83 


2.11 


0.054 


~ 1 


Fe 


(8,8) 


dxz 


-1.24 


2.01 


0.060 


~ 10- 4 






dxy 


-1.38 


2.13 


0.043 


~ io- J 



Table 3: Recapitulative table of important quantities of our Anderson models. 
Kondo temperature is so low in Fe due to the reduced broadening T and to the 
Hund coupling / that couples the two impurity orbitals. 



this prediction should in principle be amenable to experimental 
check. However, we note that our calculated Kondo tempera- 
tures are very small, which might constitute and experimental 
challenge. 

This work was sponsored under PRIN/COFIN, and also un- 
der FANAS/AFRI. We are indebted to A. Ferretti, R Lucignano, 
R. Mazzarello, and L. De Leo for early collaboration and dis- 
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Figure 1 : A schematic view of clean (4,4) and (8,8) SWNTs with their symme- 
tries. 
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Figure 4: Above: conductance as a function of energy of conduction electrons 
for Co on (4,4) SWNT, for each symmetry s - a and spin direction u - d (up- 
down); below: phase shift of conduction electrons for different symmetries s-a, 
e - o and spin directions u - d. 
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Figure 5: Predicted zero-bias anomaly for Co impurity on both (4,4) and (8,8) 
SWNTs (above) and for Fe on (8,8) SWNT (below) (g = G/G , Go = 2e 2 /h). 
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